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Abstract 

Among metal pollutants silver ions are one of the most toxic forms, and have thus been assigned to the highest toxicity 
class. Its toxicity to a wide range of microorganisms combined with its low toxicity to humans lead to the development of a 
wealth of silver-based products in many bactericidal applications accounting to more than 1000 nano-technology-based 
consumer products. Accordingly, silver is a widely distributed metal in the environment originating from its different forms 
of application as metal, salt and nanoparticle. A realistic assessment of silver nanoparticle toxicity in natural waters is, 
however, problematic and needs to be linked to experimental approaches. Here we apply metatranscriptome sequencing 
allowing for elucidating reactions of whole communities present in a water sample to stressors. We compared the toxicity of 
ionic silver and ligand-free silver nanoparticles by short term exposure on a natural community of aquatic microorganisms. 
We analyzed the effects of the treatments on metabolic pathways and species composition on the eukaryote 
metatranscriptome level in order to describe immediate molecular responses of organisms using a community approach. 
We found significant differences between the samples treated with 5 u.g/L AgN0 3 compared to the controls, but no 
significant differences in the samples treated with AgNP compared to the control samples. Statistical analysis yielded 126 
genes (KO-IDs) with significant differential expression with a false discovery rate (FDR) <0.05 between the control (KO) and 
AgN0 3 (N03) groups. A KEGG pathway enrichment analysis showed significant results with a FDR below 0.05 for pathways 
related to photosynthesis. Our study therefore supports the view that ionic silver rather than silver nanoparticles are 
responsible for silver toxicity. Nevertheless, our results highlight the strength of metatranscriptome approaches for 
assessing metal toxicity on aquatic communities. 
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Introduction 

Engineered silver nanoparticles (AgNP) are used in a wide 
variety of applications, for example as antimicrobial additives in 
textiles, as household products and in medical applications. The 
recent upward trend in production (estimated 500 t/a worldwide) 
[1] and application resulted in an increasing release of AgNP as 
well as of ionic silver into the environment as can be seen from 
elevated levels of Ag in the aquatic environment [2-5] . So far, the 
impact of AgNP, as well as of ionic silver species on aquatic 
organisms has been studied mosdy in laboratory experiments using 
single test species, sometimes even clonal cultures (e.g. Chlamydo- 
monas sp.) [6-9] . As a general trend it appears that toxicity of silver 
is due to ionic silver as the molecular toxicant [10,11]. 
Nevertheless, toxicity of AgNP is still relevant as particles represent 
a source from which Ag + can be formed continuously with 
subsequent toxic effects [7]. 

Realistic assessment of nanoparticle toxicity (mediated by their 
ionic forms) in natural waters is difficult due to the interaction of 



nanoparticles and ions with other inorganic and organic molecules 
[2]. Accordingly, it is necessary to transfer laboratory results to 
field conditions. Also, the use of single species as test organisms as 
well as analyses of single parameters such as cell numbers or 
chlorophyll content will be insufficient if community effects and 
functional diversity of ecosystems are of interest [12]. In this 
context, a metatranscriptome sequencing approach is able to 
elucidate reactions of whole communities present in a water 
sample to stressors like toxic substances [13]. Differential 
transcription of genes related to various metabolic pathways (e.g. 
photosynthesis, fatty acid biosynthesis or glycolysis) is not only 
linked to single organisms, but shows the ecological functionality of 
certain groups of taxa in a sample [14—16]. Therefore, this method 
allows detection of possible environmental hazards in a realistic 
approach, taking into account the species community as a whole. 

To the best of our knowledge, no information exists on the 
effects of silver nitrate (AgN0 3 ) as compared to AgNP on aquatic 
communities to date. Accordingly, we compared the toxicity of 
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ionic silver and AgNP by short-term exposure of a natural 
community of aquatic microorganisms in a laboratory exposure 
experiment. Since the activity of AgNP is influenced by the 
ligands, ligand-free nanoparticles are especially suitable for such 
comparisons [17]. Effects of the treatments on metabolic pathways 
and species composition were analyzed on the eukaryote 
metatranscriptome level in order to describe immediate molecular 
responses of organisms using a community approach. 

Materials and Methods 

General Experimental Set Up 

A one-day exposure experiment was conducted in June 2013 in 
a climate chamber at 16°C with homogenously distributed 
artificial day light. The intensity of the light was 60-78 |J.E 
m -2 s - ' with a 16h/8h light-dark-cycle. Approximately 150 L of 
water containing a natural plankton community from a eutrophic 
pond at the campus Essen of the University Duisburg-Essen, 
Germany, were transferred to a 200 L glass tank. The next day, 
10 L of pond water from the glass tank were filled to to nine 20 L 
plastic tanks respectively and aerated by aquarium pumps. The 
nine tanks were divided into three experimental groups (control, 
AgNOij and AgNP) with three replicate tanks each. 

Silver exposure was performed using a Ag-standard solution 
(TCP-Standard Silber, lg Ag/L, Bernd Kraft GmbH, Duisburg, 
Germany) for the AgNO^-group and a freshly laser generated 
silver nanoparticle suspension for the AgNP-group. For each 
treatment, silver was added to the water resulting at a nominal Ag 
concentration of 5 P-g/L, which was shown to be sublethal in pre- 
test experiments (see Figure SI). Monitoring of silver concentra- 
tions during exposure was performed by Ag analyses of 10 ml 
water samples taken from each tank 30 min and 24 h after the 
start of exposure; from the silver exposed groups one additional 
water sample was drawn after 5 h following exposure start. Half of 
the water samples were filtered (0.2 |J,m, cellulose acetate single use 
filter, MACHEREY-NAGEL GmbH & Co. KG, Dtiren, 
Germany) to remove organisms and organic particles. These 
samples were considered to reflect the concentration of dissolved 
silver. All water samples were acidified with 10 p.1 HNO : . 
(subboiled from 65% HN0 3 , p. a., Bernd Kraft GmbH, Duisburg, 
Germany) and were analyzed on the same day. The experiment 
was terminated after 24 h. Before exposure to silver, samples were 
taken for the metatranscriptomic sequence analysis (2.5 L of water 
containing the native plankton community) and for determination 
of water characteristics ( 1 L) from the 200 L glass tank. After 24 h 
of exposure the same sample volumes were taken from each 



treatment group. Additionally, temperature, pH, conductivity and 
0 2 -concentration were measured twice during the exposure 
period in every tank. 

Preparation of AgNP 

Silver nanoparticles were generated via laser ablation of a silver 
target in liquid aqueous medium according to Barcikowski & 
Compagnini [18] and Zeng et al. [19]. To this end, a silver foil 
(Goodfellow, 99.99%) was ablated with a Nd:YAG nanosecond 
pulsed laser (Rofin PowerLine 20E) at ^= 1064 nm with a 
repetition rate of 10 kHz and a pulse energy of 0.3 mj. The 
ablation process was carried out in a flow-through chamber with a 
volume of 1.8 ml, while the target was constantly covered with a 
liquid layer of 5.5 mm [20]. The carrier stream contained sodium 
phosphate buffer (pH 7) at an ionic strength of 50 p.M and was 
continuously pumped through the ablation chamber using a 
peristaltic pump (Ismatec ISM321C) at a constant flow rate of 
1 1.5 mL/min. To avoid the inhibition of silver ion release, which 
was shown by Grade et al. [17], the generation of silver colloid was 
carried out without additional stabilizer. Characterization of the 
sample was done by UV-Vis spectroscopy (Thermo Scientific 
Evaluation 201), recording spectra from 200-900 nm in a quartz 
cuvette (volume 3.5 ml, path length 10 mm). A significant 
plasmon resonance of silver nanoparticles (Fig. la) was detected 
via UV-Vis measurement. A narrow surface Plasmon resonance 
peak at X = 393 nm was found, which indicates the formation of 
small spherical nanoparticles, while no agglomerates (scattering in 
the NIR regime) were detected [21]. Particle size and particle size 
distributions were analyzed via analytical disc centrifugation (CPS 
Instruments Disc Centrifuge DC24000) and TEM (Philips CM12). 
TEM micrographs confirmed the formation of spherical nano- 
particles with mean particle diameters of 6 nm (distribution PDI 
0.33, Fig. lb). 

Metal Analyses and Water Chemistry 

Metal analyses of the water samples were carried out by 
electrothermal atomic absorption spectrometry (ET-AAS) using a 
Perkin-Elmer model 4 1 1 0ZL atomic absorption spectrometer 
equipped with a Zeeman effect background correction system 
(Perkin-Elmer, Massachusetts, USA). Twenty microlitres of the 
samples were injected without addition of a modifier in a pyrolytic 
graphite furnace tube with L'vov platform by the autosampler AS 
70 and run under optimized operating parameters with pyrolysis 
at 600°C and atomization at 1700°C. Calibration was performed 
by matrix adapted calibration using water of the control group 




Wavelength [nm] Diameter [nm] 



Figure 1. a) Absorption spectrum and b) size distribution data from analytical disc centrifugation of laser-generated silver 
nanoparticles with an exemplary TEM micrograph (insert). 

doi:1 0.1 371 /journal.pone.0095340.g001 
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Figure 2. Ag concentrations in the tank water (mean ± SD of three replicate tanks) as determined by ET-AAS. Initially added Ag 
concentration was 5 ug/L. Ag concentrations in the control were below detection limit. 
doi:1 0.1 371 /journal.pone.0095340.g002 



which was spiked with increasing amounts of Ag. Concentrations 
of Ag were calculated by fitting linear regression lines to the points 
defined by the spiked concentration values and the corresponding 
integrated peak areas in each sample. Correlation coefficients were 
always >0.99. 

The concentration of ammonium, dissolved phosphate, nitrate 
and nitrite were determined using Spectroquant test kits (Merck 
KGaA, Darmstadt, Germany). For the determination of the 
chloride concentration and the carbonate and total hardness 
titrimetric tests (Merck KGaA, Darmstadt, Germany) were used. 
The determination of chlorophyll a was performed according to 
DIN 38412-16:1985-12 [22]. All analyses were performed before 
and after exposure. 

Microscopical Analyses 

For the visual analysis of exposure effects, aliquots of the 
plankton samples were monitored by light microscopy. Therefore, 
protist communities were analyzed from lugol-fixed samples 
following established protocols by using Sedgewick Rafter 
chambers and Utermohl chambers [23-25]. Protists were 
analyzed at 200x magnification under an inverted microscope 
(Nikon Eclipse Ti); for small taxa 400x magnification was applied. 
Bacteria were counted from formaldehyde-preserved subsamples 
by means of epifluorescence microscopy (Nikon Eclipse 80i) after 
DAPI (4',6-diamidino-2-phenylindole) staining at lOOOx magnifi- 
cation. 

Metatranscriptome Analyses 

Following exposure, RNA was extracted from 0.2 Urn polycar- 
bonate filters using TRIzol (Life Technologies, Paisley, Scodand - 
modified). For lysis and homogenization, the cells were ground in 
liquid nitrogen in a mortar and pestie and incubated for 15 min 
with TRIzol. RNA, DNA, proteins and lipids were separated in 
phases by adding chloroform and subsequent centrifugation. The 
RNA containing aqueous phase was transferred to a clean reaction 
tube and precipitated with isopropanol. The RNA pellet was 
washed three times with 75% ethanol and afterwards resuspended 
in DEPC water. Preparation of the cDNA library as well as 



sequencing was carried out using an Illumina HiSeq platform via a 
commercial service (Eurofins MWG GmbH, Ebersberg, Ger- 
many). After quality control, one amplified short insert cDNA 
library (poly- A enriched) with an insert size of 150-400bp was 
prepared per sample, individually indexed for sequencing on 
HiSeq 2000 and sequenced using the paired-end module. In the 
following steps the preprocessing of raw reads was performed and 
the trimmed and filtered reads were subsequently mapped to the 
Uniprot database [26] for annotation. Transcript quantification 
and differential gene expression analysis was conducted thereup- 
on. 

Quality Control and Preprocessing of Sequencing Data 

The quality control tool FastQC (http://www.bioinformatics. 
babraham.ac.uk/projects/fastqc/) was used to analyze the quality 
distribution of the raw reads. Adapter sequences at the ends of the 
reads were removed using the cutadapt software [27]. Cutadapt 
was also used to trim bad quality bases with a quality score below 
20 and discard reads with a length below 30 bp after trimming. 
The amount of rRNA in the samples was determined by mapping 
the reads to the SILVA rRNA database [28] using Bowtie2 [29], a 
short read aligner that maps sequencing reads efficiently by using a 
Burrows-Wheeler transformed index. The index was built from 
the downloaded SILVA database release 1 1 1 . All quality trimmed 
reads were mapped as single-end reads against the index to 
determine the amount of remaining rRNA in each sample. Only 
the unmapped reads were used for further metatranscriptomic 
analysis. 

Mapping of Metatranscriptome Sequences 

All remaining reads were mapped to the UniProt Knowledge- 
base [26] at the amino acid level using RAPSearch2 [30]. 
RAPSearch2 uses a reduced amino acid alphabet for a very fast 
protein similarity search. We built the RAPSearch2 index from the 
downloaded UniProtKB (version May 2012) and mapped each 
single read of a read pair against the index. For each pair, the hit 
with the highest score was chosen as protein annotation. The 
mapping from UniProt IDs to KEGG Orthology IDs (KO-IDs) 



PLOS ONE | www.plosone.org 



3 



April 2014 | Volume 9 | Issue 4 | e95340 



Community Effects of Silver Nanoparticles 



o c 



I 
a. 



a. 

£ 



<J\ 0\ <J\ 



y c -z. 



W -4— 1 

1? ° 

.2 c 

0 u 

2 « 5 

*~ OJ Ul 

<u m o 



bo .id, 



E 
o ^ 

-t- -s 



[3 1] is provided by the UniProt database, and the corresponding 
KO-IDs were assigned to the reads. Mapping results were 
summarized as a count matrix of 14100 KO-IDs x 12 samples 
with the number of counts for each gene (KO-ID) in each sample. 

Normalization 

The count matrix was normalized using the "weighted trimmed 
mean of M-values" (TMM) method from the R package EdgeR 
[32]. This transforms the raw counts into counts per million 
(CPM) by normalizing for different sample sizes. 

Statistical Analysis 

To explore sources of variation in the normalized count matrix, 
we used correspondence analysis (CA), as implemented in the R 
package vegan [33]. CA is a technique that maps high- 
dimensional data onto a low-dimensional space (here two- 
dimensional plot) by singular value decomposition of the 
correspondence matrix. Each axis reveals relations between groups 
of samples and data points. Samples and data points having high 
similarity with respect to this relation have similar coordinates in 
the plot. For reasons of clarity, only the samples were depicted in 
the plots. 

For testing differential expression of genes between sample 
groups, the R package EdgeR [32] was used. EdgeR models count 
data as negative binomial distributed. The gene-wise dispersion is 
estimated by conditional maximum likelihood, and an empirical 
Bayes procedure is used to shrink the dispersion by borrowing 
information between genes. An exact test is used to test for 
differential expression between groups with a model-based 
normalization. Generalized linear model (GLM) likelihood ratio 
tests are used to determine differential expression in complex 
experiments with multiple factors. The GLM likelihood ratio test 
was applied to the metatranscriptome count data for each gene to 
account for all identified sources of variation, yielding a p-value for 
each KO-ID and treatment group. 

To visualize significant differences between groups, the p-values 
obtained from the statistical tests were plotted as a histogram. If no 
significant effect is present, the p-values follow a uniform 
distribution by definition. According to Pounds and Morris [34], 
a p-value distribution can be modeled by a beta-uniform mixture 
model, where the signal component is represented by the beta 
distribution and the null component by the uniform distribution. 
Thus a beta-uniform mixture model was fitted to the p-value 
distribution using the R package BioNet [35]. 

The significantly differentially expressed genes (KO-IDs) from 
the GLM test were used subsequently in an enrichment analysis. 
Methods from the R package iSubpathwayMiner [36] and own 
implementations were used to perform a hypergeometric test for 
each KEGG pathway. All mappings of genes to KEGG pathways 
and pathways with a significant enrichment were reported. 

Results 

Metal Analyses and Water Chemistry 

The Ag concentrations in the tank water of the different 
experimental groups were lower in the filtrated water samples 
compared to the corresponding unfiltrated samples (Fig. 2). In all 
tanks the Ag concentrations remained constant after 5 h of 
exposure, except for the filtrated AgNP samples which showed a 
slight decrease at the end of the exposure period. No differences 
were detected for physical and chemical water parameters 
between the experimental groups and their replicates. Mean ± 
SD values for each experimental group are summarized in table 1 . 
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Figure 3. Results of the morphological cell counts: A) phototrophic prokaryotes, B) heterotrophic prokaryotes, C) phototrophic 

eukaryotes, D) heterotrophic eukaryotes, E) Ceratium sp. 

doi:10.1371/journal.pone.0095340.g003 



Organismic Changes 

Microscopic analyses revealed qualitative and quantitative 
changes within the plankton community (Fig. 3; Table SI). The 
phototrophic prokaryotes (mainly consisting of Microcystis wesenber- 
gii) showed no significant changes between the two silver 
treatments (Fig.3A). The cell counts of phototrophic prokaryotes 
were slighdy higher in the silver treatments even though this was 
not significant. Nevertheless, this may indicate a similar effect of 
silver as observed for the eukaryotic phototrophs (see below). 
Discrepancies between the two silver exposure groups were 
documented in the cell counts of heterotrophic prokaryotes. The 
control and the AgNP treatment nearly had the same mean cell 
count (11.2*10 6 cells/ml and 10.8*10 6 cells/ml respectively) 
whereas the mean cell count of AgN0 3 was only 7.5*10 6 cells/ 
ml (Fig. 3B). The dominating phototrophic eukaryotes in all 
treatments were Chlorophyta, mainly consisting of the groups 



Scenedesmus / Desmodesmus and Pediastrum/ Somstrum (see Table SI). 
Both key groups showed an increase in the mean cell counts in the 
two silver treatments (Fig. 3C). One taxon which was highly 
influenced by AgN0 3 is Ceratium sp. The cell counts showed only 
13 individuals in the AgN0 3 treatment in comparison to 53 in the 
control group (Fig. 3E). 

Metatranscriptome Sequencing Results 

From the 12 environmental RNA samples, between 5.6 and 
15.2 million read pairs of 2x100 bp were obtained with good 
mean quality values of 35. After preprocessing (removal of 
sequencing adapters, low-quality parts and rRNA reads), we 
obtained 81% to 92% high-quality reads (Table 2), which were 
used for the metatranscriptome analysis. 

Each remaining read was assigned a UniProt ID and a KEGG 
Orthology ID (KO-ID). The KO-ID converts the species-specific 
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protein annotation from the UniProt database into ortholog 
groups for all proteins and functional RNAs present in the 
metatranscriptome samples, independent of the species. This 
resulted in a gene count matrix of 14100 KO-IDs x 1 2 samples, 
where each entry in the table corresponds to the number of times 
this KO-ID was identified in the sample. 

Statistical Analysis Results 

We analyzed and compared normalized counts between the 
four sample groups (raw samples "RAW", control "KO", after 
treatment with silver nanoparticles "NP", and with silver ions 
"NO3") and between the three replicate groups (suffixes "—1", " — 
2", " — 3"). Figure 4 shows a clustering of samples from a 
correspondence analysis (CA) of the normalized count matrix. On 
the first two principal axes (CA1, CA2), a clustering due to 
replicate groups is predominandy visible (Fig. 4A). The four 
different treatment groups cluster according to the second and 
third axis (CA2, CA3; Fig. 4B). 

After statistical testing for gene-wise differential expression 
between groups, accounting for the batch effect using a GLM with 
a multifactor design, we obtain p-values for each gene (KO-ID) 
and each pair of treatment versus control groups. Figure 5 depicts 
the obtained p-value distributions. For the test between control 
and AgN0 3 groups (Fig. 5A) we observed a deviation from the 
uniform p-value distribution induced by the null hypothesis. The 
deviation emerges from an enrichment of significant p-values 
oudined by the beta component in a fitted beta-uniform mixture 
model (BUM model). The beta component of a BUM model 
describes the signal in the p-value distribution, the amount of 
significant differences between the two groups, whereas the 
uniform component arises from the null hypothesis. In contrast, 
the test between control and AgNP showed no significant results 
(Fig. 5B). The p-value distribution contained no signal in the NP 
vs. control tests, as seen by the lack of a beta component in the 
BUM model. Thus, there are significant differences between the 
samples treated with AgN0 3 compared to the controls, but no 
significant differences in the samples treated with AgNP compared 
to the control samples. 

Statistical analysis yielded 126 genes (KO-IDs) with significant 
differential expression with a false discovery rate (FDR) <0.05 
between the control (KO) and AgN0 3 (NO 3) groups. The 
normalized count matrix for the significant genes is visualized as 
a bubbleplot in Figure 6, where the counts for each gene are 
displayed as filled circles. Figure 6A depicts all genes that cannot 
be annotated to a KEGG pathway, while in Figure 6B up to two 
pathways are specified for each gene. The 126 significant genes 
were used for a KEGG pathway enrichment analysis. Three 
pathways showed significant results with a FDR below 0.05. The 
enriched pathways were: photosynthesis - antenna proteins, 
carbon fixation in photosynthetic organisms and photosynthesis 
(Figure 6 and Table 3). 

Data Availability 

The raw sequence data in FASTQ_ format will be made 
available at the Short Read Archive (SRA) under NCBI accession 
number SRP040767. 

Discussion 

In the present study, we analyzed effects of silver applied as 
AgNO s and AgNP on the composition and metatranscriptome of 
a natural community of aquatic protists from a eutrophic pond in 
a 24 h laboratory exposure. To our knowledge, this is the first 
approach assessing the impact of silver on ecological functions and 
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Figure 4. Correspondence analysis of normalized count matrix. Location of samples (RAW, KO, NP, N03), each in three replicates (-1, -2, - 
3) is shown after correspondence analysis on the first two principal axes (A: axes CA1, CA2) and on the second and third principal axis (B: axes CA2, 
CA3). Roughly, the first two axes cluster according to replicate group, while the second and third axes cluster according to treatment. 
doi:1 0.1 371 /journal.pone.0095340.g004 



shifts in species composition on transcriptome level. The silver 
exposure concentration applied in this experiment was 5 Hg/L as 
either AgN0 3 or AgNP. According to pre-trials, higher concen- 
trations were found to alter the planktonic community in a way 
that was visible microscopically. For example, all eukaryotic 
organisms were dead after 24 h exposure to a silver concentration 
of 100 (Xg/L using AgN0 3 . Even at a silver concentration of 
10 M-g/L (AgN0 3 ) fewer flagellates and dinoflagellates survived a 
one-day exposure compared to unexposed controls. Thus, we 
opted for a silver concentration of 5 Hg/L to guarantee survival of 
most protist taxa, but similarly to induce effects in the community. 
This concentration is close to environmentally relevant data which 
range between 0.01 to 65 (Xg/L for freshwater ecosystems in 
Germany [37]. At all concentrations during the pre-trials AgNP 
exposure resulted in less pronounced effects than ionic silver, 
which is in line with the assumption that the toxicity and 
antimicrobial effect of AgNP is considered to be caused by 
bioavailable Ag + ions released from nanoparticles by oxidation 



[38]. Ion leaching from silver nanoparticles is a process which is 
influenced by factors such as storage conditions [39], electro- 
chemistry (influence of ions) [40] , medium additives (like proteins) 
[17], and chemicals [41,42]. Simultaneous interactions of these 
factors, which cannot be analyzed under controlled analytical 
conditions at once, provide complex kinetics of ion release. In the 
present study at least the role of biomolecules cannot be excluded 
completely. This means that not all released ions are bioavailable 
or toxic. For example, in the case of Ni-alloyed nanoparticles it 
was shown by Hahn et al. that high ion release from nanoparticles 
does not correlate with higher toxicity if the ions are bound in 
complexes by albumin or citrate molecules [43]. Further, it was 
shown that citrate deactivates silver ion release more than albumin 
[42] . Since all experiments in the present study were conducted in 
the same way and we have not registered any effects of the 
nanoparticles used, no detailed silver leaching experiments from 
Ag nanoparticles were performed. Generally, toxicity of AgN0 3 is 
higher compared to AgNP [11,3] with toxic concentrations and 



B 



"in 



CD 

Q 




'55 



CD 

Q 




0.0 0.2 0.4 0.6 0.8 1.0 



P-value 



P-value 
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inhibitory concentrations being in a range of 0.1 to 20 mg/L Ag + 
for prokaryotes. For eukaryotic cells toxic concentrations of silver 
ions range between 1 to 10 mg/L and for silver nanoparticles 



between 10 to 100 mg/L [1 1]. However, it has to be stressed that 
only little research has been conducted on toxicity of AgNPs to 
planktonic taxa [3]. Recently, Ribeiro et al. [44] described 
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Table 3. Significantly enriched KEGG pathways are shown along with their pathway ID, pathway name, the ratio of genes mapped 
to the pathway from all selected genes and the ratio of genes belonging to the pathway from all genes in all pathways, p-value for 
the enrichment and FDR. 



pathwayld pathwayName annMoleculeRatlo annBgRatio PValue FDR 

path:00196 Photosynthesis - antenna proteins 6/126 28/14039 1,49E-07 4,46E-05 

path:00710 Carbon fixation in photosynthetic organisms 5/126 35/14039 1,41E-05 1,72E-03 

path:00195 Photosynthesis 6/126 61/14039 1,72E-05 l,72E-03 



doi:1 0.1 371 /journal.pone.0095340.t003 

impaired reproduction of Daphnia magna at 1 |rg AgNP/L and 
0.5 Hg/L (Ag j after 21 d of exposure, and feeding rates being 
affected at 10 |ag/L AgNP and 2 Hg/L Ag + after 24 h exposure. 
Nevertheless, AgNP can have toxic effects that are higher than 
expected according to the concentration of dissolved ionic silver 
[7], probably due to additional effects of particles and agglomer- 
ations on cell membranes [45], depending on various factors like 
media used [46], organic molecules, light conditions and particle 
size or NP coating [9,47,48]. In the present experiment we 
detected a toxic effect on protists based on cell counts for the 
AgNO s treatments (5 Hg/L) and an effect on their gene regulation 
which is in accordance with the findings of other studies [8]. 
Effects seem to be taxon-specific, e.g. the presence of a cell wall in 
a protist species seems to buffer the negative effects of both AgNP 
and AgNO s in comparison to those protists lacking one [8,46] . In 
our study the effects were also clearly species-specific with some 
taxa such as Ceratium being affected stronger than others. However, 
on the community level we did not observe significant shifts, i.e. 
the community composition based on higher taxonomic ranks 
(orders, phyla) remained the same after one-day exposure. This 
may change following longer exposure periods as it could be 
shown in a study using a marine mesocosm that bioconcentration 
and trophic transfer of silver occurred among different taxonomic 
groups [49]. 

The AgNP treatments also showed no significant differences in 
transcriptomic response compared to the control. It may be 
assumed that only low concentrations of Ag + were available to 
cause effects within the exposure period. On the other hand, the 
same concentration of AgN0 3 showed significant differences in 
transcriptomic response compared to the control. To avoid a bias 
in the functional analysis due to the overexpression of a single 
gene, a misannotation or a false positive significance, we analyzed 
the changes between the silver treated samples and controls on 
pathway level. Pathways were reported which were significantly 
enriched in deregulated genes and therefore the effect of a single 
gene was compensated. Three pathways were significantly 
deregulated in the presence of ionic silver: photosynthesis - 
antenna proteins, carbon fixation in photosynthetic organisms and 
photosynthesis. The effects of Ag + ions on several organisms based 
on differences in transcribed genes were analyzed in previous 
studies [15] where photosynthesis was also amongst the most 
strongly affected pathways. However, genes with photosynthesis- 
related function were mostly down-regulated after exposure to 
silver in previous studies [15], but up-regulated in our study. In 
our experiment all samples were dominated by a bloom of 
Chlorophytes and no significant change of the Chlorophytes and 
heterotrophic organisms in the different treatments was observed. 
In our view, the most likely explanation for the deviating results 
concerning photosynthesis-related genes are the different concen- 
trations of silver applied in the respective studies. For example, 
Simon et al. [15] exposed Chlamydomonas reinhardtii to a AgNP 



concentration of 1 mg/L which is 200 times the concentration of 
our experiment. Silver ions are probably taken up through the cell 
membrane via the Cu (I) transporter [50], due to similar properties 
of both metals. Therefore, Cu and Ag + might cause a similar 
cellular response as the cells might not discriminate between the 
effects of these two metals. Copper ions are known to bind strongly 
to chlorophyll [51,52] and thus the efficiency of photosynthesis 
should decrease in the presence of Cu ions. To compensate such 
an effect, expression of genes relevant for photosynthesis will be 
up-regulated as we found for AgN0 3 in our experiments. Further 
evidence for such silver effects in the literature are scarce, but 
Watanabe et al. [53] reported a similarly strong binding of silver 
to chlorophyll, which may cause an up-regulation of genes related 
to photosynthesis in the presence of sublethal levels of silver. At 
higher silver concentrations direct toxic effects on metabolic 
pathways may be the dominating effects which thus decrease the 
photosynthesis rate. Such a reduction in photosynthetic yield was 
reported by Navarro et al. [7] in Chlamydomonas reinhardtii exposed 
to AgN0 3 (EC50:188nM after 1 h) and AgNP (EC50:3300M after 
1 h) and by Wang et al. [54] in Raphidocelis subcapitata exposed to 
AgNO :5 (EC 50:290 nM) and AgNP (EC50:1 1 12.63 uM) after 
4.5 h. 

Supporting Information 

Figure SI Results of the pre-test experiments. In a 

preliminary set of experiments we exposed a plankton community 
from a eutrophic pond at the campus Essen of the University 
Duisburg-Essen, Germany, to different concentrations of AgNO :i , 
as a basis for selection of silver concentrations to be used in the 
main experiment. For the pre- tests, approximately 30 mL of pond 
water were transferred into cell culture flasks and exposed to 
AgNO :5 (0ug/L, 0.01 |i,g/L, 0.1 ug/L, 1 jtg/L, 10 u.g/L, 
100 IXg/L) under the same experimental settings as described for 
the main experiment. After 24 h sub-samples were taken from the 
cell culture flasks and checked for the occurrence of living cells 
under the light microscope. Since most dead protist cells lyse 
within minutes to hours, only the living cells were counted. The 
pretest focused on the abundance of heterotrophic protists. In 
addition we used one phototrophic dinoflagellate, i.e. Ceratium sp.: 
Ceratium sp. seemed to be a sensitive indicator organism when 
exposed to silver. Further, for this species the enumeration of living 
and dead cells was possible. Therefore, the ratio of living to dead 
Ceratium sp. individuals were counted as well. Analysis of survival of 
heterotrophic protists as well as of the ratio living : dead Ceratium 
sp. showed EC 50 -values ranging between 1 to 10 Hg/L. Accord- 
ingly, we have decided to apply a Ag concentration of 5 Jlg/L. 
(DOCX) 

Table SI Microscopically cell counts of planktonic 
organisms. 

(DOCX) 
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